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We demonstrate a universal physical mechanism to probe the macroscopic quantum phase tran- 
sition based on circuit QED architecture. We found that, with certain parameters, the Josephson 
junction qubit array behaves as an antiferromagnetic Ising model in transverse field and the coupled 
transmission line resonator serves as a bosonic quantum probe. Our investigation shows that, at 
the critical point, the drastic broadening of the spectrum of the probe indicates the quantum phase 
transition. 
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Introduction- Non-analyticality in ground state en- 
ergy of a quantum many body system at critical point is 
referred to as quantum phase transition (QPT) 1], which 
is essentially caused by quantum fluctuations even at zero 
temperature. Some recent investigations have discovered 
that 1^ Q the critical behavior of a system with QPT 
can enhance quantum decoherence of its coupled external 
system. Here, the enhanced decoherence is displayed by 
a sensitive decay of the Loschmidt echo (LE) 5] at criti- 
cal point and possesses some universality in the ordered 
domain These discoveries enlighten us to propose 
a scheme to probe the intrinsic QPT phenomena of a 
system by detecting the exotic spectral structure of its 
coupled system. 

On the other hand, as a macroscopic QPT phe- 
nomenon, the superfluid-Mott insulator transition has 
been demonstrated in a macroscopic quantum system - 
the atomic Bose- Einstein condensate in an optical lattice 
0. It is natural to extend the research on QPT to some 
other macroscopic quantum systems, such as the super- 
conducting Josephson junction (JJ) array system. For 
the generic JJ array, much effort has been devoted to the 
superfluid-Mott insulator transition, see Ref.Q and ref- 
erences therein. In " qubit" regime, this system has been 
studied for some other purposes e.g., quantum 

state transfer. 

In this paper, we investigate the macroscopic QPT 
of an Ising chain in transverse field (ITF) implemented 
with JJ qubit array for the first time. We present and 
study a physical mechanism to probe its QPT with a cou- 
pled on-chip superconducting transmission line resonator 
(TLR) We find that, when the QPT occurs in the 

JJ qubit array, the spectrum of the TLR is significantly 
changed from discrete-peak structure into almost white 
noise spectrum. This drastic broadening of the spectrum 
serves as a witness of QPT. We also discuss the univer- 
sality of QPT exhibited in the spectrum structure. 

QPT model based on JJ qubit array and ID TLR - 
We consider a quantum network including N Cooper 
pair boxes (CPBs) (see Fig.©). Each CPB is a dc- 
SQUID formed by a superconducting island connected to 



two Josephson junctions. The effective Josephson tun- 
nelling energy can be tuned by the magnetic flux 
threading the dcSQUID. With proper bias voltage, the 
CPB behaves as a qubit and then JJ qubit array 
becomes an engineered "spin" chain with N 1/2 -spins. 
When the coupling capacitance Cm between two CPBs 
is much smaller than the total one Ce (e.g, in ref. [l^ . 
Cm/CY. ~ 0.05), the terms ~ o^Cm/C^) in Hamiltonian 
can be neglected and we only consider the nearest neigh- 
bor interaction in this "spin" chain. Then the JJ qubit 
array can be described by a ID ITF model with the ef- 
fective Hamiltonian 

N 
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where A = B^/B and B — e^Cm/C\ 
Coulomb interaction between nearest neighbors 
Josephson energy of CPB B^ ~ Ej cos {^x/^n) /2 with 
Ej the Josephson energy of single junction and $o = 
h/2e the flux quantum. For simplicity, all qubits are 
assumed to be identical and biased at the degenerate 
point. The quasi-spin operators cr^ ~ |0) (0| — |1) (1|, 
Cx — — 10) (1| — |1) (0| are defined in terms of the charge 
eigenstates |0) and |0) and|l) denote and 1 excess 
Cooper pair on the island respectively. A most recent ex- 
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FIG. 1: (Color on line) The schematics of our setup. A capac- 
itively coupled Josephson junction qubit array is placed in a 
ID TLR. Each qubit is coupled with the quantized magnetic 
field of the TLR. 
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periment has demonstrated the possibihty to implement 
a four-JJ-qubit Ising array jl3| . 

In our setup, as a quantum probe, a ID TLR of length 
L is placed in parallel with this JJ qubit array (see 
Fig-O) away from a distance d. Each CPB situates at 
the antinodes x = {2n+ \)L/2N (n = 0, • • • , iV - 1) of 
the magnetic field induced by the curent J in the TLR 
. Since J vanishes at the end of the TLR, the London 
equation provides the boundary condition for the electro- 
magnetic field of this on-chip resonator. Thus, the elec- 
tric field vanishes at those antinodes and the qubits are 
only coupled with the magnetic component. The mag- 



netic flux threading each dcSQUID is 



?7 (a -f I 



1 /2 

with 77 = (S/d) (hluj / L) where I is the inductance per 
unit length and S is the enclosed area of the dcSQUID. 
Here, we have assumed only a single mode of magnetic 
field with frequency w is coupled with JJ qubit array [l^ 
and a (a^) is its annhilation (creation) operator. Usually, 
Tj is small enough for the harmonic approximation [2, 
cos </)a; « 1 — 0^ and the Hamiltonian H = Hq + Hp takes 
a spin-boson form 



Hp — huja^a — g 'S^ (a^a -t- aa)) a^"-* , 



(2) 



a set of pseudo-spin operators [13 

Szk = 7I7/C +7-fc7-fc - 1, 
Sxk = i il-klk + 7-fc7i) , 



Syk = 7lfc7l-7-fc7fc- 



(4) 



They describe the pairing of quasi-particle excitations 
by 7fe. With these pseudo-spin operators, each branch 
Hamiltonian _ff„ can be rewritten as Hn = X]fe>o-^"'^^' 
where 



H'^i^ ^ Snkiszk cos2Q;„fe -I- s^k sin2a„fc) 



(5) 



with 2ank = ~ Ok, £nk = £k{K) and 8nk = dk{K)- 

Detection of QPT - We expect to detect the critical 
behavior of the J J qubit array by the coherence property 
of the TLR. To demonstrate the quantum coherence of 
a single mode electromagnetic field, a natural option is 
the correlation spectrum function S (w) = J dte~^'^^S (t), 
which is the Fourier transformation of the 1st order cor- 
relation function of the single mode field 



S{t)^{a^ (i)a(O)) 



^n|c„|^ {t)e~ 



(6) 



with the coupling coefficient g = rjEj. Here, we have 
already invoked the rotation wave approximation to ne- 
glect the high frequency terms proportional to a^a^ and 
under the condition uj >> B^, B. This approximation 
condition can be satisfied with accessable parameters in 
current experiments. For example, if we take Cs 600 
aF and Cs ~ 30 aF, L ~ 1 cm, S" 10 /im^, d ^ 1 fim 
and N = 500, then B = 1.6 GHz, Ej = 13 GHz, lu - 120 
GHz, r?~ 0.01 (IJ. 

Pseudo-spin representation for paired excitation spec- 
trum - By introducing the Jordan- Wigner transforma- 
tion ai"^ = UfXa (24*^/3 ~ 1) + 4) and cri°' = 

1 — 2cJ,Ca, Hq can be diagonalized as iJn = J^.u £ k^l^k 
by the fermionic quasi-particle operator llTlllSj 



Here, the average (■ • • ) is taken over an initial state 
l^* (0)) = iV'o) ® \G) and |i/'o) = I]„ c„ |n) is an arbitrary 
pure state of the TLR (our discussion here is also valid if 
iV'o) is an arbitrary mixed state). The decoherence fac- 
tor Dn,n-i{t) = {G \exp{iHnt) exp{—iHn-it)\G) evalu- 
ates the overlap of the wave functions under two differ- 
ent Hamiltonians i7„ and Hn-i- We also phenomeno- 
logically introduce the decaying factor exp(— F \t\) in the 
quasi- mode treatment of dissipation 0| . For strong cou- 
pling limit, g ^ T and F is about 6.3 MHz for the first 
excitation mode [lo| . 

By carrying out the evaluation of time evolution, we 
obtain explicitly the spectrum function 



S{u;) 
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with dispersion relation £/c(A) = 2B^/1 -I- — 2A cos k 
for tan0fe(A) = sink/ (A — cosfc). The ground state \G) 
of i/o describes the state without any quasi-particle ex- 
citation. 

With respect to the Fock state |n) of TLR, the Hamil- 
tonian of the whole system can be decomposed as H = 
Y^^Hn \n) {n\ where Hn = h{Xn) are defined by Eq.QJ 
with A„ = A — (2n -f l)g/B and a constant term hnuj 
has been omitted. For further convenience, we introduce 



where 
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is the sum of many Lorentzian distributions with 
the same half width at half maximum (HWHM) 

(n) 

F, but different central frequencies fl), , = 

' ^ {(ak,bk)} 

J2ki°'k£nk + bk£n~i,k)- The sum in eq.( |HIl is taken 
over all the possible configurations of combinations 
{{ak,bk) \ak,bk — i}, e.g., one possible combination is 

{(+,-)i,(+,+)2r--(-,+)^/2}. Here, = 
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rifccilCfc^^ is defined by 



(n.n—l) 



= - sin ank cos a„_ifc sin (a„_ifc - a„fc) , 

^(ri^n-l) ^ g.^^ ^^^^ ^.^^ Q;,„_ifc COS (a„-lfe - a„fc) , 
(n.n— 1) / \ 

c__|_ j, = COS a„fc COS a„_ ifc COS (Q!„_ifc - a„fej , 

c'L"_l''fc cosa„fcSina„_ifcSin (a„_ifc - Q;„fc) . (9) 

Without considering the decay of quasimodes, those 
Lorentzian hne shapes reduce to delta functions. 

The time evolution of the 1st order correlation function 
Sif) with N = 1000 is shown in Fig.© for different 
A with IV-o) = (|0> + |l>)/\/2. It can be seen that the 
decay rates for different A are almost the same except 
A = 1. This decay is induced by the dissipation of the 
quasimodes, which has the same influence for different A. 
However, near the critical point, i.e., A « 1, the decay 
is drastically enhanced. This means that there exists an 
extra strong decay mechanism related to QPT. 

To illustrate the effect of QPT more clearly, we re- 
sort to the behavior of the spectrum function S(ui). The 
numerical result by FFT is shown in Fig.© and Q. In 
Fig.®, the left panel is plotted for Ii/iq) = (|0) + |1)) /\/2 
while the right panel for j-^o) = \a) with a = 1. It can 
be seen that, for both of the two initial states, generally 
there are only one or several Lorentzian peaks centered 
at discrete frequencies in S (lo) while near the phase tran- 
sition point the spectrum of TLR gets broad and chaotic. 
As N increases, this broadened distribution at the crit- 
ical point becomes more and more smooth and tends to 
be a white noise spectrum at large N limit (see Fig.Q). 
Thus, the QPT of the JJ qubit array is featured by the in- 
tensive broadening in the TLR output spectrum. Hence, 
from the correlation spectrum of the quantum probe, we 
can infer the occurrence of QPT. 

To investigate the underlying physical mechanism for 
the behavior described above, we rewrite S (lo) as 



|i?ir^))L(c.,A(:]\r). (10) 
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FIG. 3: (Color on line) The spectrum S(iS) is shown with 
different A. The left panel is plotted for |i/>o) = (|0) + \\))l^f2 
and the right panel is for ji/io) = |a). Here A'^ = 1000. 



Here \Ef^) is an eigenvector of iJ„ with eigenvalue E\ , 



p(")=n|c„r(G|£;f))(4"-i)|G) and 
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is a Lorentzian function with the HWHM F and central 
frequency A^"'' = — £"1," In some sense, S(uj) 
measures how many different eigenvectors of i?„_i are 
needed to express one eigenvector of The more are 
necessary, the wider the support of S (uj). This observa- 
tion provides the intrinsic reason for the widening of the 
spectrum. Since g is assumed to be small perturbation 
one would generally expects that the difference between 
Hn and Hn-i is almost negligible and their eigenvectors 
are very close to each other, that is |i?|r~^^) « Si^i' 
and 
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FIG. 2; (Color on line) The 1st order correlation function 
S{t) for IVo) = (|0) + \l))/V2 is plotted with different A. 
Here N = 1000 and the time t is in the unit of 1/B. 



S{u;)^L{Lu,0,r)Y,n\cn\^ (12) 

n 

Then the support of S (lu) is very narrow and the corre- 
sponding Fourier transformation S (t) decays very slow. 

However, the above analysis is invalid at the critical 
point. Near the critical point, the property of the QPT 
system, such as ground state and long range order is sig- 
nificantly influenced by a small perturbation in either of 
the two competing terms: the Ising interaction and the 
transverse field. The seemingly very small difference be- 
tween the two Hamiltonians iJ„ and i?„_i actually has 
drastic impact on the evolutions driven by the two Hamil- 
tonians. This implies that more eigenvectors of 77„_i are 
needed to reproduce one eigenvector of i/„. Therefore, 
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FIG. 4: (Color on line) The spectrum S{ui) at the critical 
point is plotted for diflerent with IV^o) = (|0) + \1))/V2. 



more Lorentzian shapes have to be included and the sup- 
port of S{uj) becomes much broader. This in turn ac- 
celerates the decay of S{t) and acts as the extra strong 
decay mechanism related to QPT as we have noticed in 
Fig.©. 

The mechanism described in the paper can be extended 
to the case with the two-level atom as probe. The univer- 
sality similar to Ref ^ is also revealed in Fig.Q. Here, 
when A is not large, the location of peaks in S{uj) de- 
pends on both A and \ipo). However, for very large A, 
there is only one Lorentzian peak in the spectrum and 
the location of this peak is independent of A and {ipo)- 
In Fig.|(5J), we can also see that the decay envelope for 
A = 500 overlaps with A — 100. This is because the ap- 
proximation in ea. (|12|l is rigorously hold only if the J J 
qubit array is far away from the critical point. In this 
case, the spectrum exhibits universal features. 

This probe mechanism requires g S> F, which ensures 
the decoherence related to QPT is far more prominent 
than that caused by surrounding environment. But g 
also should be much smaller than the energy scale of the 
free qubit array. Otherwise, the QPT nature of the ITF 
model would be significantly changed. 

Conclusion - In this paper, with superconducting cir- 
cuit QED structure, we demonstrate a detection scheme 
for the macroscopic QPT phenomenon. By examining 
the coherent output of the coupled TLR, the quantum 
criticality of the JJ qubit array can be probed. The de- 
veloping experiments |ld ll2L [ij make our scheme to be 
potentially feasible in the near future. Concerning exper- 
imental implementation, we would like to point out only 
the case oi = is discussed here to obtain an analyti- 
cal result. But due to the unavoidable charge fluctuation 
in our system, it is hard to set the bias charge to 1/2 pre- 
cisely. Therefore, more realistic consideration reminds us 
to concern the case of Bz ^ 0, which is modeled with a 



transverse field Ising model also with a longitudinal field. 
When the longitudinal field is weak enough, this gener- 
alized model near critical point can be revealed with a 
perturbation theory. 

This work is funded by NSFC with grant Nos. 
90203018, 10474104, 60433050, and NFRPC with Nos. 
2001CB309310 and 2005CB724508. 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 
[7: 

[8 
[9 

[10: 
[11 

[12 



[is: 

[14 



[15] 

[16] 

[17] 
[18] 

[19] 



Electronic 



address: 



suncp@itp.edu. cn| 



URL: [http : / /www . itp . ac . cii/~suncp 

S. Sachdev, Quantum Phase Transition, (Cambridge 
University Press, Cambridge, 1999). 
J. Dziarmaga, A. Smerzi, W.H. Zurek and A. R. Bishop, 
Phys. Rev. Lett. 88, 167001 (2002). 

H.T. Quan, Z. Song, X. F. Liu, P. Zanardi and CP. Sun, 
Phys. Rev. Lett. 96, 140604 (2006). 

P.M. Cucchietti, S. Fernandez- Vidal, J. P. Paz, 
quant-ph/0604136 

R.A. Jalabert and H.M. Pastawski, Phys. Rev. Lett. 86, 
2490 (2001). 

M. Greiner, O. Mandel, T. Esslinger, et al, Nature 415, 
39 (2002). 

L.S. Levitov, T.P. Orlando, J.B. Majer and J.E. Mooij, 
cond-mat/0108266 

A. Romito, R. Fazio and C. Bruder, Phys. Rev. B 71 , 
100501 (R) (2005); A. Lyakhov and C. Bruder, New J. 
Phys. 7, 181 (2005). 

M. Paternstro, G.M. Palma, M.S. Kim and G. Falci, 
Phys. Rev. A 71, 042311 (2005). 

A. Walhafl, et al., Nature 431, 162 (2004); A. Blais, et 
al, Phys. Rev. A 69, 062320 (2004). 
Y. Nakamura, Yu. A. Pashkin and J. S. Tsai, Nature 398, 
786 (1999). 

Yu. A. Pashkin, et al. Nature (London), 421 , 823 (2003); 
T. Yamamoto, Y.A. Pashkin, O. Astafiev, et al. , Nature 
(London), 425, 941 (2003). 

M. Grajcar et al, Phys. Rev. Lett. 96, 047006 (2006). 
For very large A'', the energy spectrum of the cavity mode 
is quasi-continuous. In principle it is hard to single out 
one mode especially when the two systems are not exactly 
resonant. But if we takethe dissipation for the cavity into 
account, there are only some discrete Fox-Li quasimodes 
surrounded by many addtional modes. The additional 
modes induce the decay of the Fox-Li quasimode with 
decay rate P. Therefore, the single mode approximation 
is still hold. Even though, it is worth to point out that 
our proposal here is also valid for multimode field. 
Y. D. Wang, P. Zhang, D. L. Zhou, and G. P. Sun, Phys. 
Rev. B 70, 224515 (2004) 

Y.D. Wang, Z.D. Wang and CP. Sun, Phys. Rev. B 72, 
172507 (2005). 

P. Pfeuty Ann. Phys.(N.Y.) 57, 79 (1970). 

E. Lieb, T. Schultz and D. Mattis, Ann. Phys. (N. Y.) 

16, 407 (1961). 

P. W. Anderson, Phys. Rev. 112, 1900(1958). 



